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Abstract. This article provides a pedagogical review on Klein tunneling in graphene, i.e. the peculiar 
tunneling properties of two-dimensional massless Dirac electrons. We consider two simple situations in 
detail: a massless Dirac electron incident either on a potential step or on a potential barrier and use 
elementary quantum wave mechanics to obtain the transmission probability. We emphasize the connection 
to related phenomena in optics, such as the Snell-Descartes law of refraction, total internal reflection, 
Fabry- Perot resonances, negative refraction index materials (the so called meta-materials), etc. We also 
stress that Klein tunneling is not a genuine quantum tunneling effect as it does not necessarily involve 
passing through a classically forbidden region via evanescent waves. A crucial role in Klein tunneling is 
played by the conservation of (sublattice) pseudo-spin, which is discussed in detail. A major consequence 
is the absence of backscattering at normal incidence, of which we give a new shorten proof. The current 
experimental status is also thoroughly reviewed. The appendix contains the discussion of a one-dimensional 
toy model that clearly illustrates the difference in Klein tunneling between mono- and bi-layer graphene. 
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1 Introduction 



The tunnel effect of a particle going through a potential 
barrier is now a standard exercise in elementary quantum 
mechanics, which goes back to the early days of this the- 
ory |1I2| . It is usually obtained by solving the Schrodinger 
equation either approximately with the semiclassical WKB 
method or exactly for piecewise constant (square) poten- 
tials [3]. The probability for the particle to cross the po- 
tential barrier decays exponentially with the width and 
the energy height of the barrier. Thus, even if classically 
the probability to go through the barrier is equal to zero, 
quantum dynamics allows the crossing with a tiny prob- 
ability. The mechanism behind quantum tunneling is the 
possibility for a quantum particle to enter a classically for- 
bidden regions as permitted by Heisenberg's uncertainty 
principle. This is possible thanks to evanescent waves. 

However it would be considered paradoxical for a par- 
ticle to tunnel with certainty regardless of the height and 
the width of the barrier. It turns out that such an effect 
has been described theoretically by the Swedish physicist 
Oskar Klein in 1929 [4 5J for relativistic electrons using 
the three dimensional (3D) massive Dirac equation (i.e. 
the original Dirac equation describing a relativistic mas- 
sive electron [5] in ordinary space). This effect has since 
been know as the "Klein paradox". 

Recently a similar effect, though for 2D massless Dirac 
electrons, has been predicted [7_ (see also [5^ and [3]) and 
evidences of its observation in a graphene sheet were re- 
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ported |10lllll2ll3j . The latter is now known as the "Klein 
tunnel effect" or "Klein tunneling". 

The aim of the present article is to review this effect 
and to present simple derivations that can easily be re- 
produced in order to demystify it. In particular, we wish 
to emphasize that the "Klein tunnel effect" is not a tun- 
nel effect in the usual quantum mechanical sense as it 
does not crucially rely on evanescent waves and that it is 
not paradoxical (at least from the solid state physics per- 
spective). For example, there is no problem with charge 
conservation as often stated in the context of the Klein 
paradox. We also give simple physical arguments in order 
to understand the unusual tunneling behavior of massless 
Dirac particles and stress the crucial importance of the 
(sublattice) pseudo-spin conservation. 

The article is organized as follows. In section [51 we re- 
view the low-energy effective theory of valence electrons 
in graphene and then discuss properties of the 2D mass- 
less Dirac equations that are to be used in the following 
sections. In particular, we emphasize the importance of 
the pseudo-spin conservation and prove that it implies the 
absence of backscattering at normal incidence on any po- 
tential profile. Section |3] discusses the case of a massless 
Dirac electron incident on a potential step (both sharp 
and smooth) and section 2] that of a square potential bar- 
rier. In section [SI we review in detail the current status of 
experiments on Klein tunneling. Conclusions are given in 
section[B] An appendix treats Klein tunneling in ID using 
a toy model that allows us to easily compare mono- and 
bi-layer graphene. 

We mention the following reviews on the Klein para- 
dox by Calogeracos and Dombey |5I14| and more recently 
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on Klein tunneling in graphene by Beenakker [T5] and by 
Pereira et al. [T6]. 

2 Two-dimensional massless Dirac equation 
in graphene 

2.1 Low energy description of electrons in graphene 

We start by briefly reviewing the low-energy effective de- 
scription of valence electrons in a grapheire flake. Grapheire 
is a two dimensional honeycomb crystal of carbon atoms. 
The honeycomb crystal is however not a Bravais lattice. It 
is made out of a triangular Bravais lattice with a two atom 
basis (usually called A and B). This can alternately be 
seen as two triangular sublattices. As a consequence of this 
two-site basis, the electronic wavefunction is a bispinor: 
in other words, the electron carries - in addition to its 
usual spin 1/2, which we shall neglect in the following - 
a pseudo-spin 1/2 associated with its sublattice degree of 
freedom. We shall refer to it as sublattice pseudo-spin. 
The electronic band structure of graphene is usually ob- 
tained in the nearest-neighbor tight-binding model |17j. 
Conduction (or tt) electrons are allowed to jump from the 
2^2 orbital of a carbon atom to one of its three nearest 
neighbors with hopping amplitude (a.k.a. as resonance in- 
tegral) 7 = t « 3 eV. The following dispersion relation 
results (see figure [T]): 



E = (^1 + 4cos2 (kyCi) + 4 cos (fc^a) • cos {kxV^a)j 

where the lattice constant a ~ 2.46 A and the carbon- 
carbon distance is a/\/3 ~ 1.42 A. Conduction (CB, a = 
+1) and valence bands (VB, a = —1), respectively, corre- 
spond to the different signs in the above dispersion rela- 
tion (a = ±1 is called the band index); they touch at two 
inequivalent points - Dirac points or valleys K and K' - 
which are at the corners of the hexagonal Brillouin zone 
in reciprocal space (see Figure [T]). These K and K' points 
are separated by a distance ~ 1/a in reciprocal space. 

In the vicinity of the Dirac points the energy depends 
linearly on the wave number, similarly to a massless rel- 
ativistic (or ultra-relativistic) particle (see figure [T|). As a 
consequence, at low energies, and close to the K and K' 
points, the electrons can be described by a 2D massless 
Dirac eigenvalue equation: 

- ihvpS- ■ Vtp(r) = Eil^iv) (2) 

A detailed derivation of this linearized equation starting 
from the tight-binding model can be found, for example, in 
|18j . Here vp = y/i^a/ {2h) « 10® m/s is the Fermi veloc- 
ity in graphene which plays the role of an effective velocity 
of light; a = {(Tx,&y) is the 2D vector of Pauli matrices, 
■tp{r) = (-(/'^(r), -0B(r)) is the two-component (bi-spinor) 
wavefunction of the electron, E its energy and p — >■ — iftV 
is the momentum operator in the position representatioifl 

^ Note that in the low energy effective description, the 
wavevector k and the corresponding momentum p = hk are 




Fig. 1. Band structure of graphene as computed in the 
tight-binding model of Wallace [T7]. The energy E is 
shown as a function of the two dimensional Bloch wavevec- 
tor {kx, ky). In the vicinity of the Dirac points at the two 
inequivalent corners K and K' of the hexagonal Brillouin 
zoire, the dispersion relation is linear and hence locally 
equivalent to a Dirac cone (see the zoomed region). In 
undoped grapheire, the Fermi energy lies exactly at the 
two Dirac points and the Fermi surface consists of just 
two points: the valence band (VB) is filled with electrons 
and the conduction band (CB) is empty. From a band 
theory perspective, undoped graphene is therefore a zero 
gap semiconductor. However, as it is experimentally found 
that it still conducts electricity at the lowest attainable 
temperatures (with a coirductivity of the order of a few 
times the conductance quantum e^//i), it is better called a 
(zero band overlap) semi-metal. [Figure courtesy of Mark 
Goerbig]. 

The two compoirents of the wavefunction refer to the two 
atoms A and B in the unit cell. There are actually two 
such Dirac equations: one for each Dirac point or valley 
{K or K'). In the following, we only consider a single Dirac 
cone, as if we could separate the K and K' valleys. This is 
a valid approximation if intervalley scattering is unlikely, 
which is the case if the poteirtial changes are smooth on 
the lattice scale. This point is discussed in more detail be- 
low, see sectioir 1^751 The effective descriptioir iir terms of 
massless Dirac electrons is valid only for energies smaller 
than the bandwidth 7^3 eV. 

For a general reference on graphene see [TH] and for a 
pedagogical introduction see j20i2T| . 



2.2 Eigenstates of the 2D massless Dirac Hamiltonian 

In matrix notation, the two-dimensional massless Dirac 
hamiltonian is given by: 

iJkin = hvFk.a- = hvF ^ -^^ Q ^ (3) 



now defined from the K or K' point and no more from the 
center F of the Brillouin zone. The restriction to low energy 
also means that |fc| -C 1/a. 
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where hats ' denote 2x2 matrices in sublattice space 
{A, B) and the index "kin" refers to the kinetic energy. 
For simplicity, in the following, we use units such that 
ft = 1 and Wi? = 1, therefore energies and wavevectors 
are equal. As it will be useful when considering piecewise 
constant potentials, we add a constant potential to this 
hamiltonian (which simply amounts to shifting the zero 
of energy) : 

H corresponds to the total energy, iJkin to the kinetic 
energy and VqI to the potential energy, where 1 is the 
unit 2x2 matrix. From equation @ we have the eigen- 
value equation: -ffkinlV') = ^'kinlV') where fi'kin = E — Vq. 
An eigenstate 1-0) corresponds to a plane wave, which we 
write: 

^(^1^) (5) 

where k = {kx,ky) is the wavevector and the bispinor 
is given by u (respectively v), which is the complex am- 
plitude on the A (resp. B) sublattice. The corresponding 
kinetic energy is such that E'^^^ = k'^ + ky. In the follow- 
ing, we assume ky to be real and positive (i.e. a plane wave 
propagative from left y < to right y > 0). However, k^ 
can be real or purely imaginary, because fc^ can be posi- 
tive or negative. If kx is real, the wave is oscillating. If k^ 
is purely imaginary, the wave is evanescent. 



and Skin = E — Vq = ct^k"^ — , with ky — K^ > and 

a = sgn(£'— Vq). This wave decays towards increasing 
X. On the other hand, if k^ = —in 

And -Bkin ~ E ~ V{) — a\jky^ — K^,with k^ — > 
and a = sgn(i? — Vq). This wave decays towards de- 
creasing X. 

— If i?kin — (i.e. k'^ + ky = 0), again two sub-cases 
occur. On the one hand, if k^ = iky: 

On the other hand, if k^ — —iky: 

^ ^ e^'^^^e'^''^ Q (10) 

These solutions will be useful in the following when 
considering piecewise constant potentials. 

2.3 Potential steps and barriers in graphene 



2.2.1 Oscillating wave 

If kx^ > (i.e. E'^jjj > ky), the wave is oscillating and the 
wavefunction reads: 



Energy 



Ak.r 



(rife, a) 



(6) 



where a — sgn(_Ekin) is the band index (a = 1 in the 
conduction band and —1 in the valence band) and (j) is 
the angle between the wave vector fe and the x-axis such 
that tan0 = ky/k^. The kinetic energy is E^kin — E — Vq ^ 

^ + fca;^,with ky^ + kx^ > 0. The total surface A of 
the graphene sheet is taken to be one in the following 
A=l. 



2.2.2 Evanescent wave 

An interesting possibility is that fc^ be negative (i.e. E^^^ < 
ky), then kx = zLin, with k E R+. We need to consider 
two possibilities depending on Skin- 

— If -BkiiiT^O; there are two sub-cases depending on the 
sign of \kx\. On the one hand, if kx = in 



1 



(7) 




zone 2: p doped 



zone 1 : n doped 



Fig. 2. Band structure across a square potential step Vq 
(or sharp np junction). At equilibrium the chemical poten- 
tial is uniform ^{T — 0) = Ep: the Fermi level is shown 
as a dashed line. The black dot represents the electron be- 
fore and after the step: its direction of motion is indicated 
by a blue arrow. Note that its wavevector is reversed but 
not its velocity. The typical bandwidth 7 '-^ 3 eV is also 
shown. 



We consider a potential step V{r) and call w the char- 
acteristic length scale over which it varies. In graphene, 
it is possible to realize potential steps that are smooth 
{w ^ a) on the lattice scale a ^ 0.2 nm and therefore do 
not induce inter-valley scattering as the distance between 
valleys in reciprocal space is \K — K'\ ^ 1/a and the 
Fourier transform of the potential V{q) is non-zero only 
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for q <C 1/a. For such potentials, valleys are decoupled 
and electrons in graphene can be described by a single 
valley 2D massless Dirac hamiltonian H = -ffkin + V{x, y) 
|23| . In addition, as the potential varies slowly over the 
distance between two neighboring atoms {A and B), the 
potential matrix V{x, y) can be taken to be diagonal in 
the sublattice space V{x, y)l In the following, we will 
consider potentials that are translationally invariant along 
y and therefore of the form V{x)l. Therefore the complete 
hamiltonian reads: 

H = //kin + V{x)i = k& + V{x)i (11) 

In addition, the potential can be smooth {w ^ l/kp) 
or sharp {w ^ 1/kp) on the Fermi wavelength scale 1/kp — 
1/l^'kinl- If it is sharp, we can assume that it is piecewise 
constant ("square"). For example, V{x) — Vo0{x) for a 
square step (see figure^ and V{x) = Vo0{x)0{d — x) for 
a square barrier (see figure where is the Heaviside 
step function. The case of a sharp step is discussed in sec- 
tion [Xll that of a smooth step (see figure IH]) in section [23] 
and that of a sharp barrier in sectional 

Depending on the doping (the position of the Fermi 
level), the steps and barriers can correspond to different 
types of junctions. If the doping is such that the Fermi 
level lies in different bands before and after the step, the 
junction is said to be bipolar {np or pn junction). If it lies 
in the same band, the junction is said to be unipolar {nn' 
or pp' junction). Similarly barriers can correspond to bipo- 
lar (npn or pnp) or unipolar {nn'n or pp'p) junctions. In 
the following we will consider np or nn' junctions (steps) 
and npn or nn'n junctions (barriers) . 

2.4 Velocity and probability current 

For later purpose, we define the average velocity and the 
average probability current of an eigenstate \k,a). Using 
the Heisenberg equation of motion, the velocity operator 
can be evaluated as: 

v = r = -[r,H] = & (12) 

i 

From which, one can obtain the average velocity of a plane 
wave of momentum k and band index a: 

k 

V = {k,a\v\k,a) = a— (13) 
k 

Hence in the case where the band index is a = —1 - i.e. 
when the electron is in the valence band and has a negative 
kinetic energy i?kin = E—Vq < - the wavector is opposite 
to the propagation of the wave (i.e. to the velocity), which 
is quite unusual. 

The average current is obtained in the following way. 
The probability density of the state \tjj) is lip {r,t)\'^ and the 
associated (average) probability current is called j{r,t). 
The conservation of probability 

V-J--|^|V^P (14) 



allows us to define the average current as: 

j = 7/.t^^ (15) 
The average velocity is given by: 

(iplvl^p) = J (fr^/j'i&iP = J (frj (16) 

For an eigenstate |fc, a), the average current reads: 

3 = a\ (17) 

which is just equal to the average velocity as we took a 
unit system area A=l. 

2.5 Electron versus hole in the valence band 

An electron moving in the valence band (VB) should not 
be confused with a hole, which is the absence of an electron 
in an otherwise filled band (see figure |3]). Let us compare 
an electron of momentum fee in an otherwise empty VB 
with a hole that corresponds to removing of an electron of 
momentum fee in an otherwise filled VB. The electron has 
a negative kinetic energy Ef, = — /ce (as measured from 
the Dirac point), a negative charge — e, an average veloc- 
ity {vg) = —ke/ke opposite to its momentum and carries 
an electric current —e{ve). Concerning the hole, its mo- 
mentum is opposite to that of the electron fcj = —fee 0, 
its kinetic energy is positive Et = kt = —Ee, its charge is 
positive +e, its average velocity is the same as that of the 
electron (vt) — kt/kt — {v^) (therefore it moves in the 
same direction as the electron) and its electric current is 
opposite +e{vt). 

Let us now discuss Klein tunneling for an electron (ini- 
tially in the conduction band (CB)) incident on a step (see 
figure [2]). The electron is transmitted inside the step as an 
electron in the VB and not as a hole (as often stated). 
Conservation of the electric current should be enough to 
understand that point. Note that in order to keep mov- 
ing in the same direction, the electron has to reverse its 
momentum when going from the CB (outside the step) to 
the VB (inside the step). 

2.6 General conservation laws: energy E, momentum 
projection ky and ID current jx 

There are three fundamental conservation laws that we 
will keep on using when considering a massless Dirac par- 
ticle incident on a potential step or barrier that is trans- 
lationally invariant in the y direction V{x,y) — V{x)l. 

Despite the fact that on the figure[3]the hole appears to be 
at the same position (in reciprocal space) as the electron, its 
momentum is opposite kt = —fee. This is due to the fact that 
the hole corresponds to the removal of an electron. The same 
remark applies to the kinetic energy Et = — -Ee as measured 
from the Dirac point (band crossing). 
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electron in VB hole in VB 



Fig. 3. Difference between an electron and a tiole both 
in the valence band (VB). We compare an electron of mo- 
mentum fee in the VB (on the left) with a hole that cor- 
responds to removing air electron of momentum fee in an 
otherwise filled VB (on the right). The electron has a neg- 
ative kinetic energy Eg = —kg, a negative charge — e, an 
average velocity (vg) ~ —k^/ke opposite to its momen- 
tum and carries an electric current —e{ve). The hole has 
a momentum opposite to that of the electron kt = —fee, 
its kinetic energy is positive Et — kt — —Ee, its charge is 
positive +e, its average velocity is the same as that of the 
electron (vt) — kt/kt = (ve) (therefore it moves in the 
same direction as the electron) and its electric current is 
opposite +e{vt). The blue arrow indicates the direction of 
motion. 

First, there is the conservation of energy i? as a result of 
time translational invariance. Then, the momentum pro- 
jection ky is also conserved as a result of translational in- 
variance along y (i.e. parallel to the interface (s)). Finally, 
since the system is time independent, the probability con- 
servation law reads V • j = and translational invariance 
along y further implies that j{x,y) = j{x) and therefore 
that the ID current is conserved: 

jxix) — constant (18) 

2.7 Conservation of pseudo-spin and the absence of 
backscattering 

2.7.1 Pseudo-spin and chirality (or helicity) 

An electron described by the 2D massless Dirac equation 
carries a pseudo-spin 1/2 related to its freedom of belong- 
ing either to the A or to the B sublattice. In graphene, 
the electron has still an extra spin-type degree of freedom. 
It is related to its freedom of being either close to the K 
point or to the K' point in the Brillouin zone: this is called 
valley pseudo-spin. In the present paper, we do not dis- 
cuss this degree of freedom and rather concentrate on the 
sublattice pseudo-spin cr. The fact that the electron has a 
sublattice pseudo-spin is encoded in its wavefunction be- 
ing a bispinor and the hamiltonian -ffkin — k ■ & being 
a 2 X 2 matrix in sublattice space. Chirality is here de- 
fined as follows. The chirality (or helicity) operator is the 
projection of the sublattice pseudo-spin operator on the 
momentum direction: 



Its eigenvalues are C = ±1. When there is no potential 
V{r) — 0, the chirality operator commutes with the hamil- 
tonian and is therefore a conserved quantity. The hamil- 
tonian and the chirality can be diagonalized by the same 
eigenvectors: 

C\k,a) = a\k,a) (20) 

which shows that the chirality C is just the band index a 
in that case 0. 

2.7.2 Chirality factor and the absence of backscattering 

Here, we discuss an important consequence of the pseudo- 
spin, first discovered by Ando and coworkers in the context 
of carbon nanotubes [24J. Consider a massless Dirac elec- 
tron, which is incident on an impurity whose potential is 
smooth on the lattice scale such that intervalley scatter- 
ing is suppressed and the problem can be described within 
a single valley model (see section [2.3^ . The impurity po- 
tential is therefore Vimp('') ~ U{r)l For simplicity, 
though the argument can be made much more general (see 
[24J and the next paragraph), we will compute the scatter- 
ing probability using the first order Born approximation. 
It is given by 

P{e) (x\{k',a'\U{r)i\k,a)\'^ (21) 

where \k,a) and |fe', a') are the initial and the final states 
respectively and is the angle between the final and initial 
wavevectors. As the collision is elastic k' = k and a' = a. 
Therefore the only freedom in the final state is the angle 
6 = (j)k' ~ 4'k that fe' makes with fe. We are now in a 
position to compute the matrix element: 

(fe',a'|[/(r)i|fe,a) = i^^L/(fe'-fe) (22) 

where lJ{q) = j (PrU{r) ex-p{iq ■ r) is the Fourier trans- 
form of the potential U{r). Note that the transferred mo- 
mentum is q = 2fc sin(0/2). Therefore, the scattering prob- 
ability reads: 

P{0) « \U{q)\' X cos^ ^ = \Uiq)\' x i±|^ (23) 

The first term |[/(<7)p is the usual result of the Born ap- 
proximation and the second cos^ | is due to the sublattice 
pseudo-spin and is called the "chirality factor". The latter 
is just the square of the scalar product between the incom- 
ing and outgoin gbispinors: (1, e*'^'=)/V2 and (1, e*"^^' )/\/2. 
The effect of the chirality factor is quite dramatic as it kills 
backscattering (fe' = — fe): 

P{e tt) cx \U{q)\^ X cos^ | = with q ^ 2k (24) 

^ In the case where both valleys K and K' are considered, 
one finds that for an eigenstate [fc, a, ^) of the hamiltonian the 
chirality C = a x ^ is the product of the band index a and the 
valley index ^ = +1 (if if) and -1 (if iC'). 
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An intuitive explanation of this absence of backscat- 
tering is the following: if the electron tries to backscatter 
fc' = — fc it also has to reverse its sublattice pseudo-spin 
(J — !■ —<j as the pseudo-spin direction is tied to that of 
the momentum (indeed remember that away from the 
impurity i?km — k ■ &). However, the potential U{r)l 
does not act in sublattice space (it is the unit matrix) 
and can therefore not reverse the pseudo-spin. Therefore 
backscattering is impossible. This has profound physical 
consequences on the transport properties of massless Dirac 
electrons, such as weak antilocalization [25i26i27j . 



2.7.3 Conservation of pseudo-spin and the absence of 
backscattering 

We now prove that a 2D massless Dirac electron normally 
incident on a potential V{x)l can not be backscattered 
as a consequence of the conservation of its pseudo-spin 
ax- We assume that the electron is initially {t ~ 0) in 
a momentum eigenstate (kx > 0,ky = 0) and incident 
on a potential that is translationally invariant in the y 
direction {V{x, y) = V{x)). It is described by the following 
hamiltonian: 

H = kxCFx + ky&y + V{x)i (25) 

The velocity operator in the x direction is Vx — ~i[x, H] = 
ax ■ Its time evolution is given by the Heisenberg equation 
of motion 0: 

Vx = -i[^x,H] = 2azky (26) 

Here, because of the translational invariance along the y 
direction, the momentum ky is a conserved quantity: ky — 
—i[ky,H] = 0. Therefore the momentum operator along 
y does not evolve ky{t) — ky{0). If the initial state of the 
electron |'0(O)) is an eigenstate of zero momentum in the 
y direction ky{0)\tp(0)) = 0, then at any time i > 0: 

= 2(V'(0)|a,(t)A:^(0)|V(0)) =0 (27) 

which means that the velocity (or the pseudo-spin) along 
a; is a constant of the motion: 

m)\y.\m) = = +^ 

. The electron is therefore perfectly transmitted and its 
motion is exactly the same as in the absence of the po- 
tential (it is not even delayed). This shows that a (sin- 
gle valley) massless Dirac electron normally incident on 
a translationally invariant potential can not be backscat- 
tered. 

An alternative explanation for the absence of backscat- 
tering (invoking supersymmetry) is presented in [28J. 

* That the velocity operator does not commute with the 
hamiltonian is peculiar to the Dirac equation and is respon- 
sible for the so-called zitterbewegung. 



3 Potential Step 

Let us first consider a square (sharp) potential step of 
height Vq on which an electron of energy E = kp > 
is incident (see figure [2]) ■ Two zones can be defined, one 
for a: < corresponding to a kinetic energy of i?kin — E, 
another for a; > corresponding to a kinetic energy of 
iJkin = E — Vq. An analogy with an optical system will 
be made: the system is equivalent to a light beam going 
through a discontinuity between two transparent media. 
For instance, going from glass to air. Since E > in zone 
1, a = +1; whereas, in zone 2, a = sgii{E — Vq) = ±1. We 
consider an incoming electron with a given wavevector k 
in zone 1 and call k' the wavevector in zone 2. Wavevectors 
and angles are defined in figure 2] In the case where the 




Fig. 4. An electron of energy E is incident from the left 
on a square potential step of height Vq such that Vq > 
E. Angles (0 and 9t) and wavevectors in the two zones 
(before and after the step) are defined as follows: incident 
k = {kx,ky) = E{cos(j),sm(f>); reflected fe^ = {—kx,ky) = 
E{cos(tt — (j}),si'n{7: — (j))); and transmitted kf = (k'^^ky) = 
~{E — Vo){cos9t,smdt). (p is called the incidence angle. 

band index is —1 (i.e. when — Vo < 0) we saw that the 
momentum is opposite to the propagation of the wave (i.e. 
to the velocity), because the electron moves in the valence 
band. Therefore k'^ < 0, see figure 2] This gives rise to 
anomalous refraction. 

3.1 Energy below a sharp step (0 < £' < Vq) 

This corresponds to a np junction. From the conserva- 
tion of the total energy and of the momentum projec- 
tion ky, relations between angles can be deduced. The 

different wavectors are written as fe'^-') = (kx\k\^^) — 
\k''^^\{cos0j,sm9j), where E'kiiij = aj|fc^-'^| withj = {incident, 
reflected, transmitted}. Then the conserved wavevector 

projection along y reads ky — aj-Bkin.j sin^j and we ob- 
tain the table below: 



j 


incident 


reflected 


transmitted 


Band index aj 


1 


1 


-1 


Kinetic energy fi'kinj 


E>Q 


E>Q 


E-Vq<Q 


Angle 0j 


0, = 4> 


0r 


Ot 
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The equality k 



gives 



and fc^*"* — ky' gives 



(t) 



E sin ( 



-{E~Vo)sinet 



(28) 



(29) 



This last equation can be seen as an analog of the Snell- 
Descartes equation for light refraction into a medium of 
negative refraction index n [55]. To see that the refraction 
index is indeed negative, one needs to define angles from 
the perpendicular to the interface - as is usual in optics. 
For the transmitted wave ("refracted ray") in zone 2, we 
define 62 = Ot — tt and 61 = (p for the incident wave in 
zone 1. In that case both angles 9i and 62 have a modulus 
smaller than tt/2. Hence the Snell-Descartes law now reads 
nisin6'i = n2sin6'2 with ni (x E and 712 o: {E — Vq) — 
£^kin where 712 < 0. Therefore the refraction index n cx 
^-kin = ttfc is proportional to the kinetic energy and may 
be defined as n = aka ~ -Ekin/T, where 7 « 3 eV is the 
hopping amplitude (or bandwidth) and a is the lattice 
spacing. Note that the wavevector k = \k\ changes from 
zone 1 to zone 2, because the velocity vp ~ 10^ m/s is a 
constant and the kinetic energy changes from positive to 
negative. Another way to see that the refraction index is 
negative is to realize that the refraction is anomalous, or 
in other words that the refracted ray is closer in direction 
to the reflected ray than to the incident one (see Figure 
(HI), which is quite unusual. 

In optics, the refraction index n — c/v oc 1/v is in- 
versely proportional to the phase velocity v ~ where c w 
3 X 10* m/s is the light velocity in vacuum - which changes 
when going from one medium to another. As the photon 
(kinetic) energy is unchanged upon crossing the interface, 
the optical index also reads n = k/ko oc k where feg is the 
wavector in vacuum and k that in the medium. In optics, 
media with negative index of refraction - so called meta- 
materials ~ have received a lot of attention recently, see 
for example [30J. 

The next step is to write the wavefunctions in both 
zones 1 and 2 and to connect them on the interface at 
a; = 0. Here normalization of the wavefunctions is not 
needed. Using equations derived in section [51 the wave 
functions can be written in both zones as (see Figure 2]) : 



■01 



ikyU 



4'2 



= te 



1 

_|_g»(7r- 



1 



The continuity of the wavefunction is used in x 
obtain: 



1 



-r 



t and e" 



re 



te 



(30) 

(31) 
to 

(32) 



^ For the Schrodinger equation, we would have had to use 
the continuity of the wavefunction and that of its derivative as 
well. Here the two-component spinor allows the same number 
of equations just from the continuity of the wavefunction. 



which gives: 



pi4> 



-i(p 



and t 



+ e-^'P 

p — i4> pif^t 



(33) 



We now wish to obtain the transmission probability 
T from the amplitudes r and t. This requires using the 
conservation of the ID current, see equation (fT8)) . which 
here reads: 



jx [incident 

O cos (h — 



■ jx [reflected] — jx [transmitted] 



rj^ cos( 



= -Wcosl 



and therefore: 



l = |r|' 



, cos 9t 



(34) 



(35) 



This is the probability conservation law 1 = R + T. It 
allows one to identify the transmission T and reflection R 
probabilities as: 



T 



cos Of 
cosdi 



t|2 and R 



(36) 



and we finally obtain: 



T 



cos cp cos 



cos' 



( <P+Bt 



and R 



(^) 



/ <P+Ot 



where the transmitted angle is: 



(37) 



TT = arcsm 



sin ( sin <h] 

\Vq-E V 



(38) 



These equations are equivalent to the Fresnel formulae in 
optics. Note that cos < so that T > as it should. 
The transmission probability is plotted in Figure [S] 

The transmission coefficient vanishes beyond a certain 
critical angle <j)c, which is is defined as sin ipc — {Vo — E)/E 
(see figures [5] and Beyond this angle, an evanescent 
wave is created in zone 2 and a total reflection is observed. 
A critical angle only occurs when Vq/2 < E, otherwise 
the conditions to have an evanescent wave tt^-^ sin (t> > 1 

E „ 



and 



Va-E ' 



< 



Vo-E ' 

1 are never satisfied. This is similar 



to the optical phenomenon of total internal reflection. For 
instance if a light beam goes from a glass-like medium 
with an index ni to an air-like medium of index n2 < ni 
then using the Snell-Descartes law nisin(?i = n2siii92, 



® The reflection coefficient R is always equal to jrj^, whereas 
the transmission coefficient T is not necessarily given by \t\'^. 
In the present case of the potential step T = / 
There is a mistake precisely on that point in ^29) . As an il- 
lustration of the fact that T is generally not given by jtj^, we 
consider an extreme example where the transmission proba- 
bility vanishes T = although t 7^ 0. If in zone 2, the ki- 
netic energy vanishes i?kin — E — Vo = 0, there is an evanes- 
cent wave when (j) ^ 0. Then using Eq. © and fTS}, we can 
show that r = exp{i24>) and therefore T = 1 — |r|^ — 0, while 
t = l + exp(i20) / 0. 
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Fig. 5. Polar plot of the transmission probability T{(j)) 
for the potential step (values of T are shown on all axis) 
for several diniensionless energies e = E/Vo between and 
1 and with the incident angle cf) running from -tt/2 to tt/2: 
e = (blue), e = 0.25 (purple), e = 0.5 (beige), e = 0.75 
(green) and e — 0.95 (black). Note that for e = 0.75 there 
is a critical angle above which T is strictly equal to zero. 



there exists a refracted beam only if the incident angle 9i 
is smaller than a critical angle 6c = arcsin (n2/ni). 

If the electron arrives at normal incidence (0 = 0) the 
probability to go through is 1: 



r(<^ = o) = i 



(39) 



This is due to the "absence of backscattering" discussed 
previously ( see section [2.7.21 and Ref. [24j) and is a con- 
sequence of the pseudo-spin conservation of the massless 
Dirac electron. This (at first) surprising result is most of- 
ten referred to as the "Klein tunnel effect" although it is 
not a genuine tunnel effect in the quantum mechanical 
sense. Indeed it involves no classically forbidden region 
and no evanescent wave. It is a consequence of (1) the 
existence of negative kinetic energy states (valence band) 
in the step, that match the positive kinetic energy states 
(conduction band) outside the step and (2) of the conser- 
vation of pseudo-spin which permits the transition. This 
second point is quite important. The chirality factor acts 
as a selection rule. Indeed the square of the scalar prod- 
uct between the incoming and the transmitted bispinors is 
[1 -I- aa' cos{9t — (f>)]/2 = 1 as a = +1 = —a', dt = tt when 
(j) = 0. However, it is not because states are available in 
the step at the matching energy that the transition will 
necessarily happen. The case of the bilayer graphene (with 
its massive chiral electrons) is illuminating in this respect 
[7]. In a graphene bilayer, the carriers are also described by 
a bispinorial wavefunction, but with a Berry phase of 27r 
(instead of tt for the monolayer) . The low energy effective 
hamiltonian is 



bilayer 



2m- 



(40) 



for a single valley, where to* is the effective mass of the 
carriers. As the band structure is also that of a gapless 
semiconductor - although with parabolic bands 

E = ±{hkf/{2m*) 

- there are also states of negative kinetic energy avail- 
able in the step. Here, however, the pseudo-spin conserva- 
tion forbids the interband transition at normal incidence 
T((/) = 0) = 0. Indeed the chirality factor Q in this case is 
[1 + aa' cos 2(6*4 - 0)1/2 = as a = +1 = -a', 6*4 = tt 
when = 0. The ID case is studied in an appendix. 

Some special set of energy E and incidence angle (j) are 
worth mentioning. For instance, li E = Vq/2 one has (see 
also HI): 

T = cos^ = 1- {ky/kpf (41) 



and \i E (because Ot — tt): 



T = 



2 cos(j) 
1 + cos (p 



(42) 



These two cases are quite interesting (see figure [S]). There 
is no critical angle, and the transmission is always quite 
large unless (j> becomes really close to ±7r/2. There is a 
slight preference for normal incidence but no true colli- 
mation effect. This will later be compared to the case of 
a smooth step. 



3.2 Energy above a sharp step (Vb < E) 

This corresponds to a Tin' junction. For a; < the kinetic 
energy is E and for x > the kinetic energy is E — Vq > 
0. Here we will consider both Vq > and Vq < 0. The 
wavefunctions are: 



Jkyy 



J 1 

V+e 



+ e 



-ikxX 



1 



1 



V'2 = te 

and the Snell-Descartes relation reads: 
Esin(t)={E- t/o)sin6't 



(43) 
(44) 

(45) 



This leads to a new system of equations that can be solved 
to give: 

r — - 



t = 



and one eventually obtains 



jCOSf^t COS0COSf/t 



cos< 



,s^ (^) 



(46) 
(47) 

(48) 



^ Note the factor of two difference in the cosine. This is a 
consequence of the Berry phase being 2ti rather than tt. 
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where: 

Ot = arcsin ( — — — sin </) ) (49) 
\ tj — / 

When Vq > 0, this case corresponds to refraction indices 
ni > and n2 > (equivalent to usual optical materials: 
the refraction is normal) with ni > ?i2, so that there is 
also a critical angle given by sin^c — {E — Vb)/£'. This 
is typically similar to a light beam going from glass to 
air. In classical mechanics, there is a similar phenomenon. 
Consider a non-relativistic particle of mass m and energy 
E = p^/(2m) incident on a potential step of energy Vo 
such that Vo < E. The conservation of energy /{2m) = 
p'^/(2m) + Vq and that of parallel momentum p'y — py 
imply that p'\l(2m) = p1./{2m) — Vq. The particle will 
therefore be reflected if p' ^ < because a negative kinetic 
energy is classically forbidden. This defines a critical angle 
of incidence (/)c given by sin (pc = \/{E — Vo)/E, where the 
angle is defined by tan = py/px- We note that the non- 
relativistic refraction index n cx \/£kin \p\ differs from 
the ultra-relativistic one of graphene n cx -Ekin = Oi\k\. 

When Vo < 0, this case corresponds to refraction in- 
dices ni > and n2 > with ni < n2, so that there is 
no critical angle. This is typically similar to a light beam 
going from air to glass. 

Using the Schrodinger equation (parabolic band ap- 
proximation) , the fact that a particle with an energy above 
a potential step (_E > > Vb) and with an incident ve- 
locity in the same direction as the step force {F = — VF) 
can be reflected is known as a "quantum reflection", be- 
cause it can not be understood from a classical perspec- 
tive. It is particularly spectacular for perpendicular inci- 
dence (0 = 0) and low energy E ^ \Vq\ where the trans- 
mission probability vanishes T — > 0. In the present case 
of a massless Dirac electron, quantum reflection is absent 
as T((/) = 0) = 1 for any energy E > Q. This is again a 
consequence of the absence of backscattering. 



3.3 Smooth potential step 

We now turn to the case of a smooth potential step, which 
was recently considered by Cheianov and Falko in the con- 
text of a np junction in graphene [5]. Take a smooth po- 
tential step of height Vq, which occurs over a distance 2w 
(Figure ini). Smoothness roughly means \p w (below, 
we discuss a refined criterion) . The potential is taken to be 
V{x) = sign(a;)Vb/2 when |a:| > w and V{x) — Fx when 
|x| < w, with F = Vq/ (2w). We consider a symetric situa- 
tion: that of an incoming massless Dirac electron with an 
energy which is half that of the step. This energy is here 
E — Ep = Oasa result of the shift in the zero of energy - 
indeed V{x < —w) = — Vb/2. The electron has a wavevec- 
tor {kx,ky) — fci?(cos (/), sin 0) where kp — Vb/2. A nor- 
mally incident electron fcy = is perfectly transmitted as a 
consequence of the absence of backscattering. We focus on 
the case ky ^{) which gives rise to a classically forbidden 
zone close to V{x) = E (i.e. x = 0). It can be found from 

the conservation of energy ^/kx{xy + k'^ + Fx = E = 




Vo 



Fig. 6. Band structure across a smooth np junction. The 
potential height is Vq and its width is 211). The energy of 
the incoming electron (Fermi energy) is half that of the po- 
tential step. The zero of energy is chosen to coincide with 
the Fermi level. Arrows indicate the direction of motion 
before (blue) and after (pink) the step. 

and the requirement that kx{x)'^ < 0. This defines the 
region |a;| < Ic = ky/F. The size of the classically forbid- 
den zone is 2lc — 2w sin cj). An electron with ky ^ that 
wishes to go through the step has to tunnel via evanes- 
cent waves to the other side across the classically forbid- 
den zone. In the semiclassical approximation, the prob- 
ability amplitude to tunnel is given hy A fa e^^ where 
S — J^l" dx kx (x) is the (reduced) action, which is imag- 
inary when kx{x) = i|fca;(2;)| ~ ^\J^'y ~ (Fx)'^. The tun- 
neling probability is then T — \A\'^. The typical kx{x) is 
kx{0) — iky, therefore S ^ iky x 2lc and T ~ g-4/c„ic 

More precisely, S — ikylc ^_^du\/\ — v? = iTrkylc/2 and 
the probability is [5]: 

T{(j)) W e^^^v^': — e^^^l/P — g-7rfcj=.tusin^ ^gQ^ 

This result is valid for a smooth step and for incidence 
angles (j) not too close to ti/2. It satisfies r(0) = 1 and 
T{(f)) 1 for incidence angles such that |0| <C 0o = 
F/nk'^p = 1/ ^/nkpw - where (po is called the collimation 
angle -, and then rapidly goes to zero for oblique incidence 

> (j>Q. This shows that a single np junction has a colli- 
mation effect as it focuses the electronic fiow by allowing 
the transmission of only the trajectories that are close to 
normal incidence |29) . This result should be compared to 
T — cos^ (j) found previously in the case of a sharp step 
(see figure[7]). Apart from small incidence angles, chirality 
seems not to play an important role here. This tunneling 
is similar to interband tunneling in a Zener diode, see e.g. 
[31j. In a semiconductor, one can have a junction between 
an electron-doped and a hole-doped region (a so-called np 
junction). If the voltage difference across the tunnel junc- 
tion is strong enough, it is possible for electrons to tunnel 
from the conduction to the valence band across the junc- 
tion. The only difference here is that the semiconductor is 
2D, gapless and the bands are linearly dispersing. 

The transition between the abrupt {w — )■ 0) and smooth 
{w — >■ cxi) steps limit is treated in [32J. Actually, a step is 
considered smooth when the semiclassical approximation 
is valid. This happens when the typical tunneling action 
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-1.5 -1.0 -0.5 
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Fig. 7. Transmission probability T((/i) across a poten- 
tial step for an energy which half that of the barrier. T 
(on the y axis) is plotted as a function of the incidence 
angle </> (on the x axis) between — 7r/2 and 7r/2 for a 
sharp (dashed blue, T = cos^ (^) and a smooth step (red, 
T = ^-^^Fws\^^ ,i> ^-^^^ ^ ;^Q)^ rpj^g colhmation effect 
is clearly seen in the smooth step case, with a character- 
istic angle 0o = '^l\J'!^kpw ~ 0.2. 



is large \S\ ^ ka:{0)lc ~ kywsinip — kpuisiii ip ^ 1. The 
smoothness criterion is therefore kpiu sin^ (/) ^ 1. It not 
only depends on kpw but also on the incidence angle 0. 
Close to normal incidence, any step becomes sharp (see 
figure lb in Ref. [32]). If kpw < 1, the step is sharp at 
any incidence angle. If kpw > 1, the step is sharp close to 
normal incidence and smooth close to grazing incidence. 
In the "sharp step angular region" (|(/)| ^ l/^/irkpiv), 
the transmission is high and corresponds to the absence 
of backscattering; whereas in the "smooth step angular 
region" > l/y^irkpuu), transmission is exponentially 
suppressed and occurs via evanescent wave inter-band tun- 
neling. The signature of chirality for the smooth step is 
therefore in the colhmation effect: only electrons suffi- 
ciently close to normal incidence are transmitted across 
the junction. 

When the energy of the incoming electron is not half 
that of the barrier, the transmission probability becomes 
llQi: 



T{(f>) sa exp(-7r- 



2k p-^ 



CF2 



k 



Fl 



(51) 



where kpi (resp. kp2) is the Fermi wavevector in the zone 
1 (resp. zone 2). The symetric case considered above cor- 
responds to kpi = kp2 = kp. 



4 Potential barrier 

Let us now consider a square potential barrier [7] , see also 
|33| . We distinguish three zones: zone 1 for x < 0, where 
the potential is equal to 0; zone 2 for < x < rf, where 
the potential is equal to Vq; and zone 3 for x > d, where 
the potential is again 0. And study two cases: < E < Vq 
(coresponding to a npn junction) and E > Vq > (nn'n 
junction). 



zone 3: n doped 




Vo 



zone 2: p doped , 

d 



^ ZONEl 


ZONE 2 




Z0NE3 
























x=0 




eB 


x=d 




kinetic Energy : E 


E-V„ 




E 





Fig. 8. Top: Band structure across a sharp npn junc- 
tion. An electron of energy < i? < Vq is incident 
from the left on a square potential barrier of height 
Vq and width d. Note the position of the Dirac cone 
tips in the three zones. Bottom: definition of the an- 
gles and the wavevectors in the three zones: incident 
k — {kx,ky) — £'(cos sin (/)), reflected kr = {—kx,ky) — 
E{cos (tt — (j)), sin (vr — <f>)), transmitted inside the barrier 
kA = {k'x,ky) = — (£^— Vo)(cos sin (?^), reflected inside 
the barrier ks = {—k'^,ky) = —{E — Vn) {cos 9b, sin 9b) 
and transmitted kt = k. 

4.1 Energy below the barrier (0 < i? < Vq) 

This corresponds to a npn junction in graphene. Using the 
equations of section [5] the wave functions can be written 
in the three zones (see figure [5]): 



^1 



■02 



Ae 



1 

1 

pi8A 



re 



Be 



-ik^x 



1 



■03 = igikyy^ik^x 



1 

gi0 



(52) 

(53) 
(54) 



With the continuity of the spinors in x = and x = d, 
the following system is obtained: 



1 



-r + A + B 



2# ^ r^'^<t> _ ^g^flA 



Be' 



= -te'''-''^ 



^e^fe^d + Be 
= -te^l'e'''-'^ - Ae^'^^e'^'-'^ + Be'"'^' 



-, — ik'd 



(55) 
(56) 
(57) 



e-if^'.d (58) 



Pierre E. AUain, Jean-Noel Fuchs: Klein tunneling in graphene 



11 



where A, B, r and t are complex amplitudes to be de- 
termined. The easiest method to solve this system is by 
substitution and only r is worth computing 



= -2e*'^sin(fc^d) X 



(59) 



+ 9 a) + e*''-'' cos (0 - 9 a) + 2i sin {k'^d) 
The transmission coefficient follows from T = 1 — IrP: 



cos^ (f) cos^ 9 a 



cos^ (^cos^ 9a cos^ {k'r^d) + siv? (fc^c?)[l -I- sin6'^ smcj)]^ 

(60) 

^ cos^ (b with the dimen- 



where k'^d = —2-kI^\ — 2e + e"^ cot 
sionless barrier width I = Vod/{2TrhvF) — Vod/{2Tr) and 
the dimensionless energy e = E/Vq. Naturally, from the 
conservation of ky we can obtained the Snell-Descartes 
law: 

Esm(l) = -{E -Vo)sm9A (61) 
Several cases are worth investigating: 

4.1.1 Low energy (E < Vq) 

The incoming electron has an tiny energy compared to Vq 
and the transmission coefficient takes a simpler form |7|: 



T 



1 — cos^ (k'^d) sin (p 



(62) 



4.1.4 Critical angle (1/2 < e < 1) 

As with the potential step, once a critical angle defined 
as sin(/)c = (Vq — E)/E is reached, an evanescent wave 
is present in the zone 2. But contrary to what happened 
with the step, there is no total reflection. Indeed, the exis- 
tence of a second discontinuity between the zone 2 and the 
zone 3 allows the wave to be transmitted through the bar- 
rier with a reduced amplitude, just like in the Schrodinger 
tunnel effect: passing through is classically forbidden but 
quantum mechanically permitted via an evanescent wave. 
This is therefore a genuine tunnel effect. The condition to 
have an evanescent wave is that the wavevector projec- 
tion along X in zone 2 is purely imaginary: /c^^ < 00. As 
k'^ d^ — (27r^)^(l — 2e-f cos 0), the condition to have an 
evanescent wave is: : 



2e- 



cos ( 



< 



(64) 



There are no evanescent waves for an electron with an 
energy below Vb/2. The critical angle is plotted in figure 
[5] in a kind of "phase diagram". 




4.1.2 Grazing energy [E ^ Vq ) 

Here, the electron arrives exactly with the energy of the 
barrier. The transmission is entirely via evanescent waves 
(except exactly at = 0). When E^kin = 0, fc^ = iky and 
therefore 1/ cos9a — jii'kinl/fcx — ^-'^d tan6'yi — ky/k'^ — 
—i. Therefore equation (|60l) becomes: 



oscillating waves 



COS 



cosh (kyd) 



sm 



(63) 



where we used that 

cos^ (ikyd) — cosh^ (kyd) and sin^ (ikyd) = — sinh^ (kyd) 

. Note that it is not obvious a priori that a formula ob- 
tained for oscillating waves remains valid in a regime of 
evanescent waves. 



Fig. 9. "Phase diagram" for the wave nature inside the 
barrier plotted in the {(f>, e) plane. On the y-axis is the 
dimensionless energy e = E/Vq, on the x-axis is the in- 
cident angle running from -7r/2 to 7r/2. The white area 
is the zone of oscillating waves (classically allowed re- 
gion; possibility of resonances), the blue area is the zone 
of evanescent waves (classically forbidden; possibility of 
true tunnel effect via an evanescent wave). The two re- 
gions are separated by the critical angle line (pc such that 
sin 0c = (1 - e)/e. 



4.1.3 Normal incidence (0 = 0) 

If the incident angle is zero (the angle is taken from the 
a;-axis), the transmission coefficient is exactly equal to 1 
regardless of the length d and the (energy) height Vq of 
the barrier. This perfect transmission at normal incidence 
is again due to the conservation of the pseudo-spin lead- 
ing to the absence of backscattering (see section [^77)) . In 
particular, it is not an interference effect between the two 
interfaces at a; = and x = d. For such an effect, see 
section H.l.Sl on Fabry-Perot resonances. 



4.1.5 Fabry-Perot resonances (0 ^ 0) 

We now restrict to oblique incidence (0 7^ 0) and consider 
multiple interferences effects. A potential barrier can be 

* Note that the condition for having an evanescent wave in 
the barrier is the same as that for total reflection on the po- 
tential step and defines the critical angle. As runs between 
— 7r/2 to 7r/2, the function sin increases with and cos > 0. 
Therefore > 0c with sin 0c = i^<t4>sin0 > sin 0c = ^-j^, 
hence £^ sin^ > (1 - e)^ = 1 + - 2e^l - 2e + cos < 
0, which is precisely the same as 
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seen as a double interface (at a; = and x — d) and as 
the analog of a well-known optical system: a Fabry-Perot 
interferometer. The cavity is the region inside the barrier, 
which can accommodate oscillating waves - especially at 
£ < 1/2, see Fig. [HI Accordingly, the incoming wave might 
interfere with itself between the two interfaces (at x = 
and X = d) in zone 2. If the waves interfere constructively 
transmission resonances will occur where r((/) 7^ 0) = 1. 
The condition of such resonances - also known as tunnel- 
ing resonances - is [7J: 

k'^d = TT X integer (65) 

which is just the condition that a half-integer (N/2) num- 
ber of wavelengths (27r/fc^) along x fits in the cavity (i.e. 
inside the barrier) of size d: (integer/2) x (27r/fc^) = d. As 
k'^d = —Ivl^/y — 2e + e'^ cos^ the resonance condition 
involves the energy e, the length of the barrier / and the 
angle (f> and reads: 

2l\/\ — 2e + e^ cos-^ (j) = integer (66) 

It defines specific angles 0„ ^ such that T{4)n) = 1. 
These Fabry-Perot resonances are responsible for the petal- 
like shape of T(0) when plotted as a function of the inci- 
dence angle (p at fixed energy e and barrier width I (see 
figures [inland El and Ref. L2j)- 
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Fig. 11. Plot of the transmission probability T as a func- 
tion of the incident angle and the dimensionless energy 
e = E/Vq for a fixed dimensionless width of the barrier 
1 = 1. The energy e varies between and 10 and the angle 
(f) between — 7r/2 and +Tr/2. Fabry-Perot resonances are 
clearly visible. 




Fig. 10. Polar plot of the transmission coefficient T{ip) 
for the potential barrier at fixed energy e and width I. 
The two sets of parameters are the same as in [7] namely: 
(e = 0.41519, I = 4.85, blue) and (e = 0.291038,'/ = 6.91, 
red) . The petal structure is clearly seen and correspond to 
Fabry-Perot resonances. Note that T(0) = 1 as a conse- 
quence of pseudo-spin conservation and independently of 
£ and I. The black line indicates unit transmission. 



4.2 Energy above the barrier {E > Vq > 0) 

We now focus on the situation where an electron is in- 
cident with an energy larger than the one of the barrier 
(corresponding to a nn'n junction) . After proceeding with 
the same kind of computation as for the other case, we ob- 
tain a transmission coefficient which has the same expres- 
sion as (|60l) and (pT|) with the replacement 9a —Oa- 
Within this energy range, most of the properties of the 
transmission coefficient are retrieved. In the high energy 
limit (i? 3> Vb > 0) limit, we find the transmission prob- 
ability: 



T = 



cos 



1 — cos^ k'd sin^ < 



(67) 



4.3 Grazing energy (E' = Vo > 0) 

Both limits E — > Vq^ have the same common value, there- 
fore T{E = Vo) = T{E = Vo"^): 



T : 



cosh {kyd) — sin 



This result plays a key role in the calculation of the two- 
terminal conductance of a ballistic undoped graphene sheet 
using the Landauer formula (see [M] and [23). The reser- 
voirs are assumed to be made of doped graphene and the 
sample is undoped. In order to ensure that the reservoirs 
conduct much better than the undoped graphene sheet, 
one takes the limit of a very large doping of the reservoirs 
(many transverse channels in the reservoirs): this is equiv- 
alent to sine/) — > 0. In this limit (0 — 0), the transmission 



Pierre E. AUain, Jean-Noel Fuchs: Klein tunneling in graphene 



13 



coefficient becomes: 



5.2 Ballistic versus diffusive regime 



T: 



cosh (kyd) 



(68) 



This corresponds to a genuine tunnel effect through the 
sample (except a.t (p — 0). Indeed, when kyd ^ 1, T « 
4exp(— 2fcyC?) <C 1, which is of the expected semiclas- 
sical form with a tunneling probability proportional to 
exp(— 2S'i), where Si — kyd is the imaginary time action 
for the classical path inside the barrier. 

To conclude this section on the potential barrier, we 
mention that the case of a smooth (trapezoidal) barrier 
has been treated in [36]. See in particular the discussion 
of Fabry-Perot resonances. 



5 Experiments 

Below we describe the current status of transport exper- 
iments designed at observing the Klein tunnel effect in 



The samples are usually small: the distance L between two 
measuring contacts is typical less than 1 /im. This is done 
in order to be as close as possible to the ballistic regime, in 
which collisions on impurities can be neglected (note that 
our complete discussion of Klein tunneling assumed that 
we could neglect the effect of disorder). To know whether 
this is indeed the case implies to compare the mean free 
path /,„ with the step size 2ii;, the barrier width d and the 
sample size L. Depending on the situation, it is possible 
to have a sample which is, for example, globally diffusive 
{Im ^ L), that the barrier is diffusive (/„ ^ d) but that 
each step can be described as being ballistic (/„ ^ w). 
The mean free path in a good graphene sample is typi- 
cally on the order of /,„ ~ 100 nm. The effect of disorder 
on a graphene np junction is discussed in These au- 
thors find that the transition between ballistic {(3^1) 
and diffusive (/3 ^ 1) regimes is controlled by a single di- 

mensionless parameter /3 = Idn/dxl/n^^"^ , where dn/dx is 



graphene and rely mainly on the following papers fl0'll'12'13f|ie density gradient right at the junction and : 



but mention also other relevant experiments [37 38 39 40J . 
For reading convenience, in this section, we restore the 
units of h and v p . 



5.1 Backgate and topgate 

Graphene samples are usually equiped with a backgate 
that allows one to electrically control the doping of the 
sheet through an electric field effect (similar to a capac- 
itor) [51]. In exfoliated samples, the backgate is usually 
made of doped silicon and separated from the graphene 
sheet by ~ 300 nm of silicon dioxide (dielectric) . Electro- 
static potential steps and barriers can be made by using 
an additional topgate. The simultaneous use of a backgate 
and a topgate allows one to control independently the 
energy of incoming electrons (Ep) and the step/barrier 
height (Vb). The distance between the topgate and the 
graphene sheet roughly gives the step size 2w. Steps made 
in this way are smooth on the lattice scale {w ^ a, as 
an example 2w ~ 80 nm [Tn]) and until recently also 
smooth on the Fermi wavelength scale {w > 1/kp, typ- 
ically 1/kp ~ 10 nm). Sharp steps {w < l/kp) should 
soon become available thanks to rapid progress in fabri- 
catior{^. The typical barrier width d is in between 100 nm 
and 1 /xm. 

Transport measurement are performed on these sys- 
tems, in which the two-terminal resistance is measured as 
a function of the backgate and topgate voltages. Roughly 
speaking, the backgate controls the Fermi energy and the 
topgate sets the barrier properties (height Vb, width d and 
step size 2w). One may in addition apply a magnetic field 
in order to measure the magneto-resistance. 

^ Reaching small kp values is limited by the presence of 
so-called electron-hole puddles close to the Dirac point in a 
graphene sheet. These inhomogeneities mean that locally the 
Fermi wavevector is never really zero but has a finite min- 
imal value . The cleanest suspended samples correspond to 
l/kp ~ 100 nm. 



- where /i™ is the mobility - roughly gives the density of 
impurities. This prediction was confirmed experimentally 
112J. 

Low temperature is also needed in order to have co- 
herent propagation of the electrons (no decoherence) . This 
corresponds to our assumption of treating the electron as 
ideal matter waves, rather than classical particles. The 
coherence length in a graphene sample is typically > 
1 /im when the temperature is below 4 K. 



5.3 Smooth junctions and poor screening 

Potential steps realized up to now were generally of the 
smooth type. However, the slope F of the potential step 
V{x) « Fx right at the junction (when V{x) = Ep) is 
not properly estimated as Vo/2w (as we did in section 
13. 3p . Such an estimate relies on assuming perfect screen- 
ing in the graphene sheet, which is not correct. Indeed, 
close to the bipolar junction, screening in graphene is very 
poor because it corresponds to the crossing of the Dirac 
point where the density of states vanishes. Linear Thomas- 
Fermi screening would predict no screening at all, as the 
inverse screening radius vanishes. It is however possible to 
study non-linear screening in this regime. This was done 
by Zhang and Fogler for a ballistic np junction [IT. They 
find that the slope F is strongly enhanced compared to the 
naive estimate Vo/2w (typically by a factor of 10). This 
effect is important to take into account when comparing 
theory and experiment [12113] . 



5.4 Evidences for the observation of Klein tunneling 

5.4.1 Resistance of a smooth ballistic np junction 

The first experiments designed at observing Klein tunnel- 
ing in graphene all measured the resistance across a npn 
junction as a function of topgate and backgate voltages 
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|10I11I12] . Because of the presence of disorder, the resis- 
tance of the barrier was found to be correctly described as 
the sum of the resistance of two smooth np junctions (in 
series). Each np junction was in the ballistic regime but 
not the whole barrier. As the momentum of an electron 
is not conserved during its motion between the two inter- 
faces, cavity type resonances are not possible, see e.g. [44J. 
We therefore consider these experiments as testing the re- 
sistance R across a smooth step {np junction) rather than 
a barrier. 

In order to set the stage, we give typical values of rel- 
evant quantities. The barrier height Vb ~ 0.1 eV, width 
d ~ 300 nm and step size 2w ~ 100 nm; the Fermi wave- 
length l/kp ~ 10 nm such that kpw ~ 5 (smooth step); 
the mobility ^ 1000 — 10000 cm^/V.s, the mean free 
path ^ 30 — 100 nm and the length between measuring 
contacts L ^ 1.3 — 5 /im. 

These measurements probe the average transmission 
across the junction. Indeed, the two terminal conductance 
G = 1 / R is given by the Landauer formula 

ch. 

where the sum is over transverse channels labelled by ky = 
kp sin (p and with transmission probability T{ky), W is the 
sample width (not to be confused with the step size w) 
and the factor 4 accounts for valley and spin degeneracy 
in graphene. 

As an example, we consider a symmetric np junction: 
the Fermi energy is at half the potential step. On the one 
hand for a smooth step of potential slope F at the junc- 
tion, the transmission is T{ky) — q~^^'"p'^v/^ (see section 
[51^ and therefore |8j: 

On the other hand, for a sharp symmetric step, the trans- 
mission is T = cos^ 0=1 — [ky/kpY (see section [Q]) and 
therefore the conductance is: 

^sharp ^ kpW 

Zh TT ^ ' 

For comparison, when there is no step, the transmission 
is perfect T{ky) = 1 and the conductance is 

Gno stop — 4:-, (72) 

h n 

As G"4fP = 2Gno step/3 > G™°°*'', we conclude that 
the sharp np junction is almost transparent, while the 
smooth np junction is highly resistive and only lets the 
electrons close to normal incidence through (coUimation 
effect). Note that from the measurement of G™°°*'' and 
Gno step it is possible to estimate the collimation angle as 

G™°°"VGno stop ~ 0o/2. 



The resistance measured across smooth np junctions 
as a function of the topgate voltage (controlling the bar- 
rier properties) was found jllll2| in agreement with the- 
ory provided the ballistic regime is reached (/3 ^ 1) [12] 
and the poor screening in graphene close to the junction 
is acounted for [33]. Quite counter- intuitivelj0, the resis- 
tance across a smooth junction was measured to exceed 
that predicted in a purely diffusive model (i.e. excluding 
chirality effects) but to agree with the ballistic prediction. 
The signature of Klein tunneling was precisely found in 
this collimation effect that only allows electrons close to 
normal incidence to go through the smooth junction. 

Although the average transmission was found to agree 
with the prediction of Klein tunneling across a smooth 
np junction, the angular dependence of the transmission 
probability T(0) was not seen. In particular, from a mea- 
surement of the average transmission, it is not possible to 
tell that perfect tunneling occurs at normal incidence or to 
precisely measure the collimation angle. A measurement 
of the resistance across a sharp ballistic junction would 
be closer to revealing perfect tunneling as its resistance 
is predicted to be only 3/2 times larger than in the ab- 
sence of a step. In conclusion, these first experiments gave 
indirect evidences of Klein tunneling in graphene bipolar 
junctions. 



5.4.2 Conductance oscillations and magneto-resistance 
across a ballistic npn junction 

A second type of experiment was performed in order to 
test Klein tunneling more directly, which relied on quan- 
tum interferences between two pn interfaces. A graphene 
npn junction was realized with a narrow topgate such that 
the whole barrier was in the ballistic regime ^13j. Indeed, 
the mean free path was estimated as Im ^ 100 nm larger 
than the barrier width d < 100 nm. The step of size 
2w ~ 30 nm was smooth compared to both the lattice 
spacing and the Fermi wavelength 1/kp ^ A nm. The mo- 
bility was Urn ^ 5000 cm^/V.s, the distance between mea- 
suring contacts L ^ 3 /im and the typical barrier height 
Vq ~ 0.3 eV. Two main observations were made on this 
system. 

First, oscillations in the conductance as a function of 
the top gate voltage revealed that the whole npn junction 
was (for the first time) in the ballistic regime. These oscil- 
lations where interpreted as interferences due to multiple 
reflections between the two pn interfaces (see the section 
14.1.51 on Fabry-Perot resonances) . The topgate allows one 

A widespread misconception about Klein tunneling is that 
it should systematically allow electrons to go through any bar- 
rier with a high probability. If this was true, it would indeed 
seem counter-intuitive that an experimental evidence for Klein 
tunneling comes from measuring an increase of the resistance. 
Note that this increase is here defined with respect to a dif- 
fusive model and not to a situation in which Klein tunneling 
would be simply turned off. And that the evidence for Klein 
tunneling in such an experiment on smooth junctions is in the 
collimation effect. 
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to tune these transmission resonances and to span the in- 
terference fringes. Note that, due to perfect tunneling at 
normal incidence, Fabry-Perot resonances are only possi- 
ble for oblique (0 7^ 0) trajectories. 

Second, applying a perpendicular magnetic field, resis- 
tance measurements revealed a half-period shift in these 
Fabry-Perot fringes above a critical magnetic field ^ 0.3 T. 
The interpretation is as follows [48J. Reflectionless trans- 
mission at normal incidence (0 = 0)- in other words, 
Klein tunneling - also means that the reflection ampli- 
tude r undergoes a tt phase jump when the incident angle 
tp goes from positive to negative valu^Hl- At zero magnetic 
field, two consecutive (non-normal) reflections on the two 
pn interfaces occur with opposite angles and 02 — —'Pi- 
A weak magnetic field bends the electronic trajectories. 
Above a critical field, trajectory bending becomes suffi- 
cient to make the two consecutive reflections occur with 
the same incident angle 0i — 4)2- This suddenly adds tt to 
the phase accumulated by an electron between two reflec- 
tions and shifts the interference fringes by half a period. 
The observation of this half-period shift (see figure 3 in 
|13j l is therefore a direct evidence of perfect tunneling at 
normal incidence. 

6 Conclusion 

In conclusion, we compare Klein tunneling to the standard 
tunnel effect and outline what are the crucial ingredients 
necessary for its occurence. We also show that the phrase 
"Klein tunneling" has different meanings. Eventually we 
give indications for further reading. 

The standard tunnel effect across a barrier is an in- 
traband transition of a Schrodinger electron via evanes- 
cent waves across the classically forbidden zone (the bar- 
rier). It gives a tunneling probability at normal incidence 
which is roughly T(0 = 0) ~ g-2Kd .^^j^gj-g typi- 
cal wavevector ( in the x direction, perpendicular to the 
barrier) inside the barrier and d is the barrier width. The 
probability therefore decays exponentially with the width 
and the energy height of the barrier (through k). 

Klein tunneling is the name given to the interband 
transition (say from a conduction to a valence band) across 
a step or barrier of a massless Dirac electron. It relies (i) on 
having negative kinetic energy states available in the step 
or barrier (matching energy of the incoming electron) and 
also (ii) on the pseudo-spin conservation that may allow 
or not the transition (roughly speaking this is a kind of 
selection rule given by the overlap of the bispinors outside 
and inside the barrier; its modulus square is the so-called 
chirality factor). 

In section 13.11 we showed that, in the case of a sharp 
potential step, the reflection amplitude is r(0) = (e**^ + 
e''*')/(e''^ - e"'**) where Ot = -k + arcsin(£; sin (^/(Vo - E)). 
In the limit were </> = ±77 with 77 — >■ C*", it follows that r(±r;) f» 
e±™/2r;Vo/(2(V'o - E)). Therefore, there is a tt phase jump 
(in the reflection amplitude) when the incident angle changes 
sign: Argr(»7) — Argr(— 77) ~ tt when — >■ 0"*". More gener- 
ally, for arbitrary r), one can show that Argr(77) — Argr(— 77) = 
■K + 2\q + arcsin(£;sin77/(Vo - E))] = tt -I- 0(77). 



Depending on the precise situation Klein tunneling 
may refer to different physical situations and mechanisms. 
Here we distinguish four situations encountered in the 
present article (by default, it is usually the first case which 
is meant): 

1) At normal incidence on a potential step that is 
translationally invariant (along y), there is perfect inter- 
band transmission without evanescent waves: 

r(0 = o) = i 

This is a consequence of the absence of backscattering 
due to pseudo-spin conservation. It is not a genuine 
quantum tunnel effect. See section [2.7.31 

2) At oblique incidence on a sharp step, there can be 
T < 1 interband transmission without evanescent waves. 
For example, for an electron incident with an energy which 
is half that of the step: 

r(0 7^ 0) = cos^ 4) 

The chirality factor is playing a role here. See section [01 

3) At oblique incidence on a smooth step, there is in- 
terband tunneling (via evanescent waves). For example, 
for an electron incident with an energy which is half that 
of the step: 

T{(j) 7^ 0) « g—^'iVFkp sin^ 4>/F 

where F is the potential gradient at the bipolar junction. 
This is a genuine quantum tunnel effect when 101 > 0o = 
y/ F/ {nhvpkp) and chirality only plays a role close to 
normal incidence leading to collimation of the electrons 
(|0| <C 0o)- See section 

4) At oblique incidence on a square barrier, there can 
be Fabry-Perot resonances. These are transmission res- 
onances due to the multiple interferences of oscillating 
waves between the two interfaces: 

T(0„ 7^ 0) = 1 when k'^d = 7rn 

See section 14.1.5 1 

We end this article by giving a list of references for 
further reading on topics not covered in the present re- 
view. Klein tunneling has also been studied theoretically 
for massive Dirac electrons in "gapped graphene" |3S], for 
massive chiral electrons in a graphene bilayer with |46j or 
without [7J band gap and in a deformed honeycomb lat- 
tice [37] . It was also investigated in monolayer graphene in 
the presence of a magnetic field |8I48I49| or a superlattice 
[50 51]. The motion of a wavepacket through a barrier in 
graphene is discussed in [16]. 
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also M. Biittiker for encouragements and S. Gueron, B. Huard 
and P. Carmier for useful comments on the manuscript. This 
work was realized during the internship (June- July 2008) of 
P.E.A. at LPS Orsay in partial fulfillment of his master (Ml) 
degree at the Universite Paris-Sud, France. 



16 



Pierre E. AUain, Jean-Noel Fuchs: Klein tunneling in graphene 



A Klein tunneling in one-dimension 

As an illustration of the importance of pseudo-spin con- 
servation, we consider Klein tunneling in ID and study 
two toy-model hamiltonians. First consider the following 
massless Dirac hamiltonian - which we call the "ID mono- 
layer": 

+ V{x)l (73) 
with hvp = ^ where vp is the Fermi velocity. The pseudo- 
spin is conserved as [ax,Hm] = 0, and thus also the veloc- 
ity operator Vx = —i[x,Hm] = ax is a conserved quantity 
O If the electron is initially in a velocity eigenstate (say 
such that Vx = +1), then {vx{t)) — +1 at any t > 0. 
Therefore the motion of the electron in the presence of 
the potential is exactly the same as in its absence (the 
motion is not even delayed) . This is a strong consequence 
of the absence of backscattering. 

To understand the physical meaning of the pseudo- 
spin, we consider the case when the potential is absent. 
An eigenstate of the hamiltonian is: 

where kx is the momentum, -Ekin — sigii(E]^in)\kx\ is the 
energy and = ±1 is the eigenvalue of the pseudo-spin 
ax - The latter can also be written as ax = sign(A;a;)sign(i?kin) 
or ax = E]iin/kx and is therefore the direction of motion 
(+1 for right movers and -1 for left movers). If an elec- 
tron is initially a right mover, it will remain so even if it 
encounters regions of arbitrary non-zero potential. 

In the presence of the potential, it is also possible to 
find the eigenstates of the hamiltonian ([751) , see . Per- 
forming a unitary transformation the hamiltonian can be 
written as Hm — kxa^ + V{x)l. The eigenvalue equation 
decouples in two equations ^idip±/ dx = [E — T^(a;)]V'±, 
which are easily solved to give the following eigenvectors 

V.^,.Jx) ^ -LQJe^-T'^-'I^-^K)] (75) 

at any eigenenergy E. Each energy level is doubly degen- 
erate as ax = ±1- This clearly shows that, for a scalar 
potential V{x)\ at any energy, the eigenstates are delo- 
calized as |-0£;.cr^ (a;)p = constant. 

Next, consider the case of massive chiral electrons (a 
kind of "ID bilayer" toy-model) described by the following 
hamiltonian: 

Hb = klax + V{x)l (76) 

with h^/{2m*) = 1 where m* is the effective mass (see 
equation (|40p for the corresponding 2D hamiltonian). The 
pseudo-spin is again a conserved quantity. However, the 
velocity operator is not, as it is now given by Vx = -'i[x, H^] - 
2kxax and [vx,Hb] = -2iaxdV/dx ^ 0. When V{x) = 0, 
the eigenvectors of Hb are: 

= e-^- (77) 

Among Dirac equations, this is peculiar to the ID case, 
that features no zitterbewegung. 




/ \ 
/ \ 

' * X 

Fig. 12. Klein tunneling in one dimension. An electron in- 
cident from the left on a sharp potential step (the blue ar- 
row indicates its direction of motion) . Top: the "ID mono- 
layer" case, in which the pseudo-spin corresponds to the 
direction of motion: the branch of right movers is in red 
ax = +1 and that of left movers in dashed green ax — — 1- 
Bottom: the "ID bilayer" case, in which the pseudo-spin 
corresponds to the band index: the conduction band is 
in red ax ^ +1 and the valence band in dashed green 
ax = —1. In both cases, conservation of pseudo-spin im- 
poses the direction of motion after the step (indicated by 
a pink arrow). 



with the corresponding eigenenergies i?kin = axk"^, which 
shows that here ax — sign(i?kin) can also be seen as the 
band index. When V{x) is non zero, we can perform a uni- 
tary transformation to rewrite the hamiltonian as Hb — 
k'^&z + V{x)l. The eigenvalue equation then decouples 
in two ID Schrodinger equations ^d'^ip±/dx'^ — [E — 
V{x)]ip±. As for a generic potential, all states of the ID 
Schrodinger equation are localized [53J , it follows that the 
eigenstates of the "ID bilayer" hamiltonian are also local- 
ized. This is the opposite conclusion to the "ID monolayer" 
case. Here the conservation of pseudo-spin leads to local- 
ization. 

The conclusion that we draw on Klein tunneling across 
a step is twofold: the transition in the step is possible if 
(1) there are states available in the step at a matching 
energy (i.e. states of negative kinetic energy) and (2) if 
the pseudo-spin conservation permits such an inter-band 
transition. The latter provides a kind of selection rule, re- 
'flecting whether the appropriate matrix element for the 
inter-band transition vanishes or not. This matrix ele- 
ment (squared) is usually called the chirality factor and is 
given by the overlap of the incoming bispinor {l,ax)/ \/2 
and the transmitted bispinor (1, cr^)/-\/2, which in ID is 
(1 -I- axa'x)/2 = 6cr^,a'^- In the case of the monolayer, 
the pseudo-spin is the direction of motion ax = -Ekin/fcx 
and therefore the transition occurs with unit probability. 
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Whereas in the bilayer case, the pseudo-spin is the band 
index — sign(£'kin), which would obviously change in 
an inter-band transition a'^ — —ax, which is therefore 
strictly forbidden. These results are reminiscent of the 
2D case at normal incidence where T = 1 for the mono- 
layer and T = for the bilayer [7], and the correspond- 
ing inter-band chirality factors are (1 — cos7r)/2 = 1 and 
(1 — cos(27r))/2 — respectively, see section [XT] 

Note that at oblique incidence, the 2D massless case 
is quite unlike its ID counterpart. In fact, for a potential 
with y translational invariance, there is an mapping be- 
tween the massless 2D case at oblique incidence {ky ^ 0) 
and the ID case of Dirac electrons with a finite mass. 
Indeed, using the conservation of ky, the 2D eigenvalue 
equation of a massless Dirac electron [—idxdx — i^ydy + 
V{x)l)il}{x,y) = Eip{x,y) becomes a ID equation for a 
massive Dirac electron: 

{—ia-xdx + m&y + V{x)l)ip{x) = Eip{x) (78) 

where ijj{x,y) — ip{x) exp{ikyy) defines the ID wavefunc- 
tion ip{x) and m = ky is the mass 0. This ID equa- 
tion is actually that originally considered by Klein with 
V{x) = VoO{x) [4J and by Sauter with V{x) = Fx \Mi- 
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